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Abstract 

We give a derivation of the dispersion law e(p) — h?p^ /2m + V{p) — 
V(0), where V{p) is the Fourier transform of the pair interaction poten- 
tial V{r). (The interaction between particles Xi and Xj is V^(|a;i — x-j\).) 
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1 INTRODUCTION 

In the famous paper [T] , Bogolyubov obtained the spectrum (energy levels) of a 
nonideal Bose-Einstein gas. 

The formula for the spectrum reads 



Xp^-h'vp + e{p), e{p)^^(^^ + V{p)^ -V^ip). (1) 

where v is the fluid (gas) velocity and V{p) is the Fourier transform of the pair 
interaction potential T^da;^ — Xj\) for particles i and j. 



For small p, the term e{p) behaves as jply y(0), which corresponds to sound 
(phonons), provided that V{0) > 0. The last condition can be rewritten as 
Io° ^('') means that the repulsion for small r exceeds the attraction 

occurring in helium for large r. The Fourier transform V{p) tends to zero 
as p — > oo, and the leading term is of the order of p'^/2. The point po of 
minimum of the radicand, where the p-derivative vanishes, is called the roton 
part of the spectrum s{p). 
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In the same journal issue where Bogolyubov's paper was pubhshed, Landau 
gave a general argument based on experiments and showing that the curve e{p) 
should have exactly this general form, i.e., be linear near p = 0, have one point 
of maximum and then a point of minimum, and then tend to infinity. 

This curve is called the Landau curve (the 1962 Nobel prize). 

The superfluidity of helium-4, which is a Bose gas, was discovered by Kapitsa 
in 1930 (the 1978 Nobel prize). 

However, while the superfluid state in the Bose case is caused by short-range 
repulsion, it is known that the superfluid state has to be caused by long-range 
attraction, which enables the formation of Cooper pairs (the 1972 Nobel prize). 

The superfluidity of a Fermi fluid (helium-3) was discovered experimentally 
in the early 1970s by the American scientists David Lee, Douglas Osheroff, and 
Robert Richardson (the 1996 Nobel prize). 

Landau wrote a year after publishing his result: "Recently, Bogolyubov has 
succeeded in determining the energy spectrum of a Bose-Einstein gas with weak 
interaction between particles in the general case by an ingenious application of 
second quantization" [H p. 43]. 

Specifically, Bogolyubov used Dirac's idea: since the number of particles is 
large, one can deem the commutator of the creation and annihilation operators 
in the Bose statistics to be small. The so-called Bogolyubov asymptotic (u, v)- 
transform leads to the same result. 

In the Fermi statistics, one considers the anticommutator of the creation- 
annihilation operators, so that the above argument fails to provide passage 
to "c-numbers." This passage is only possible for pairs of fermions. 

The theory of ultrasecond quantization developed by the present author 
permits one to obtain an even simpler formula, based on long-range attraction, 
for a Fermi gas. Namely, 



A„ = —h pv 



2^2 



2m 



V{p) - ViO) 



(2) 



This formula is the subject of the present paper. 

The Landau curve also arises in this case, and the main role is played by the 
long-range attraction provided by the interaction potential. To show this, we 
first study the behavior of V{p). 

Remark . The second-quantized Schrodinger equation identically coincides 
with the iV-particle Schrodinger equation up to a unitary transformation. We 
mean that the second-quantized equation is a convenient representation of the A'^- 
dimensional system of Schrodinger equations. This representation (like the 
p-representation, the interaction representation, etc.) is "identical" up to a 
unitary transformation. It contains no approximations and is, in this sense, 
an identity transformation. In just the same way, ultrasecond quantization, 
introduced by the author, coincides up to a unitary transformation with the N- 
particle Schrodinger equation for symmetric and antisymmetric solutions. For 
some reason, this new identity scares off theoretical physicists so much that they 
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begin to cross themselves, saying "keep away from me," and publicly burn the 
author's books as well as his likeness. 

However, Bogolyubov's asymptotic (u, u)-transform, which is valid for bosons, 
coincides in the fermionic case with the ultrasecond quantization method only 
for Bardeen's interaction potential. This transformation is approximate rather 
than exact and gives a wrong pair interaction asymptotics. Apparently, it is 
only after a trap experiment with a Fermi gas, similar to the one carried out in 
the remarkable paper [3] , is conducted for the ultrasecond quantization method 
that theoretical physicists will believe in the method. 

2 "GENERAL POSITION" THEOREM 

In mathematics, the notion of "general position" for points of a surface, say, 
with respect to projection onto a coordinate plane, is well known and was in 
particular used in an essential way by Arnold, Varchenko, and Gusein-Zade. In 
physics, this special case of general position means in particular that the optical 
focus of a lens refracting light rays so that they all converge at a same point is 
not a generic situation — once the lens is moved slightly, the focus spreads into 
a small light spot formed by a network of intersecting caustics. 

On the plane, a perpendicular drawn to the abscissa axis is not in general 
position — after a small rotation of the coordinate axes, this line is no longer 
projected into a single point on the abscissa axis. If the first two derivatives are 
zero at a point of a smooth curve, then the curve is not in general position with 
respect to a small rotation of the axes. One can rotate the axes in such a way 
that, in a majority of cases, the smooth curve will not have such points. Then 
one speaks of a curve in "general position." 

As to a point where the derivative is zero, a small rotation destroys this prop- 
erty. Then one speaks of a "generic" point of a smooth curve. The derivative at 
a generic point of a smooth curve is nonzero. In particular, the derivative of a 
generic smooth curve on the half-line is nonzero at the origin. If this is not the 
case, then a small rotation of the coordinate axes destroys this property. Thus, 
a small rotation does not affect the property of the derivative to be nonzero at 
the origin. The following elementary theorem holds. 

Theorem . Let the original pair interaction potential U{x — y),x,yizM.^, be ra- 
dially symmetric and equal to V{r), where r = i^i ~ Vi)^ andlhiir^o V{r)r ^ 0. 

Suppose that the three-dimensional Fourier transform V{\p\), \p\ > 0, is in gen- 
eral position at the origin \p\ =0. Then 

lim p'^V{p) ^ 0, lim r^V{r) ^ 0. (3) 

p — ^oo r — *oo 

Proof . Since V{^) is spherically symmetric, it follows that so is the three- 
dimensional Fourier transform V{^). 
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By passing to the spherical coordinates, we obtain the relation 

Vir)^- r ViOsm{^r)^dC (4) 

^ J — oo 

Since F(|^|) tends to infinity at the rate of const/|^p, it follows that 

, const 

V{r) = + 0(1) 

r 

as r — > 0. The converse is also true. 

The condition of general position implies that 

V'iO) ^ 0. 

One can readily see that 

y'(0) = --5- lim rV(r). (5) 

47r ^ — *c)o 

Indeed, by replacing the integral (|4]) by the integral over the finite inter- 
val [—A, A], where A is sufficiently large, by making the change of variables 

= 77, by integrating by parts, and by passing to the limit as r — > oo 
and A — s- oo, we arrive at the assertion of the theorem. □ 

In fact, the interaction potential is not exactly spherically symmetric in view 
of the complicated structure of the molecule. In principle, one should take into 
account the influence of all particles forming the molecule. If, nevertheless, we 
treat the molecule as a whole and approximate the interaction potential by a 
spherically symmetric function, we should at least assume that the potential is 
"generic." 

Corollary . One has 

V{p) - V{0) - -Att\p\ lim r^V{r). (6) 

r — ^oo 

as IpI — !■ 0. Since the potential V{r) is attractive as r —i- oo, it follows that 

- 4:77 lim r'^Vir) = c> 0. (7) 



r — >oo 



As p ^ OO, the expression under the modulus sign in ^ tends to infinity at 
the rate of h^p^/2m, since V{p) ^ oo at the rate of l/|pp. 

Thus, we see that formula ([2]) gives the Landau curve for a Fermi fluid. 

Remark . Consider the Lennard-Jones potential and add a small smooth ra- 
dially symmetric potential to it. The probability that the new potential still 
decays at the rate of l/r^ is zero with respect to the entire set of potentials. 
With probability 1, the decay will be at the rate of There is no known 

physical law that would prohibit this, and the existence of such a law is unlikely. 
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3 ULTRASECOND QUANTIZATION 



Ultrasecond quantization for problems in quantum mechanics and statistical 
physics was introduced in [4]-[7]. Recall the notation and main facts for the 
case of quantization with respect to "particle-number" pairs and with respect to 
pairs of particles. Quantization with regard to pairs will permit us to take into 
account pair correlations of particles when constructing the asymptotics. The 
ultrasecond quantization space is the bosonic Fock space ^ [5]. Let 6"'"(x,s) 
be the creation operator for particles with number s, and let h~(x,s) be the 
annihilation operator for particles with number s in =^ [S]. Next, let B^{x, x') be 
the creation operator for a pair of particles, and let B~ (x, x') be the annihilation 
operator for a pair of particles in this space. These operator satisfy the fermionic 
anticommutation relations 

{b-{x,s),h+{x\s')} ^5,s'5{x- x'), {S±(x,s),&±(x',s')} = 0, 
[B-{xi,X2)B+{x\,x'^)] = 5{xi - x'^)5{x2 - 4), xa), 4)] = 0, 

= [h^(x,s),B^(x\,x'2)]=Q. 

(8) 

For bosons, 6 and 6^ satisfy similar commutation relations. 

Next, $0 is the vacuum vector in ^ with the following properties: 

r(x,s)$o = 0, B-{xi,X2)^o=0. (9) 

The variable x ranges over the three-dimensional torus Li x L2 x L^, which will 
be denoted by T. The variable s = 0, 1, . . . is discrete; it is called the number, 
or the statistical spin. Each vector $ G can be uniquely represented in the 
form 

* = EE ^ • • • 51 / • • ■ / . . . dxfe dyi . . . dy2M 

fc=OJ\/=0 ■ '81=0 Sfc=0'^ 

X $fc,M(a;i, Si; • • • ; a;fe, Sk] yi,y2; ■ ■ • ; y2M-i,y2M) 

X si) • • • 6+(xfc, Sfc)i3+(yi, 2/2) • • • B+ {y2M-i,y2M)'i>o, (10) 

where the function $fc,Af (xi, si; . . . ; Xfc, s^; yi, 2/2; • ■ • ; 2/2Af-i, 2/2m) is symmetric 
with respect to transpositions of the pairs {xj, Sj) and {xi, Si) as well as of the 
pairs {y2j-i,y2j) and (2/21-1, 2/21)- In the bosonic case, one introduces the sub- 
space ^t^M™ of vectors $ such that ^k',M' = for (fc', M') ^ {k, M) and $/c,m 
is a symmetric function of the variables xi,X2, . . . ,Xk,yi,y2, ■ ■ ■ , y2M- In the 
fermionic case, one in a similar way introduces the subspace -^^W^™ 
tors $ such that ^k'.M' = for (fc', M') ^ {k, M) and ^k,M is an antisymmetric 
function of the variables Xi,X2, . ■ ■ ,Xk,yi,y2, ■ ■ ■ , y2M- The orthogonal projec- 
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tion onto ^I^m"^ in J? has the form [l]-[7] 



fifeS""" = TWn ^---^ I ■■■ I d^i---dxkdyx... dy2M 



si=0 Sfc=0' 

X 6+(a::i,si) • • • 6+(a;fe, Sfe)_B+(j/i, 1/2) ■ ■ ■ {y2M-i,y2M) 
(b {xi,si)---b {xk,Sk) 

X B^{yi,y2) ■ ■ ■ B^{y2M-i,y2M)) 



exp ^— J dxh^(x,s)b (x,s) — jj 



dydy' B+{y,y')B-{y,y') 
(11) 



where Symm^^ ^^y_^ y^^^ is the operator of symmetrization with respect to the 
variables Xi, . . . , Xk, j/i, . . . , y2M and the operators b'^{x, s), b~{x, s), B~^{y, y'), 
and B~ {y, y') are Wick order [H] . The orthogonal projection onto in 
has the form [?] 



nt," = ^---^ ... dxi...dxkdyi... dy2M 

' si=0 Sk=0-^ 

X 6+(a;i, Si) • • • b^{xk,Sk)B^{yi,y2) ■ ■ ■ -B+(?/2A/~i, 2/2m) 
X Asymm^^_^^y^ ^^^^(5"(a;i,si) ■■■b-{xk,Sk) 
X B~{yi,y2) ■ ■ ■ B~{y2M-i,y2M)) 

^^^p(-£/ dxb+{x,s)b-{x,s)- JJ dydy' B+ {y,y')B~ {y,y')^ , 

where Asymm^^ ^^y_^ y^^^ is the operator of antisymmetrization with respect 
to the variables xi, . . . , Xk, j/i, . . . , y2M- From now on, the operators b'^{x, s), 
b~{x,s), B^{y,y'), and B~{y,y') are Wick ordered unless specified otherwise. 

In what follows, we consider a system of N identical particles on the torus T. 
We assume that the Hamiltonian of N bosons or fermions has the form 

,2 W N N 

= E ^- + E E " ^0- (12) 

j=i j=i i=j+i 

According to [4J, the ultrasecond-quantized Hamiltonian corresponding to this 
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operator in the bosonic case has the form 

Hb^J2J2 / ■ ■ • / dxi...dxkdyi... dy2M 

fc=OA/=0 ■ ■ si=0 Sk=f) 

X Si) • • • 5+(xfe, Sk)B^{yi,y2) ■ ■ ■ {y2M-i,y2M)Hk+2M 

X Symni^^_,^j^^...y^^^ (h~ {xi, Si) ■ ■ -b' {xk, Sk) 

X B^{yi,y2) ■ ■ ■ B^{y2M-i,y2M)) 

xexpf-^ / dxb+{x,s)b^{x,s)- // dydy B+{y,y')B^{y,y') 

(13) 

In the fermionic case, the corresponding operator Hp is expressed by a similar 
formula with Symm replaced by Asymm. By analogy with (|12p and (|13p. we 

assign an ultrasecond-quantized operator A 0] to every iV-particle operator 

1 1 
-r ( 2 2 d d 

An xi, . . . ,XN]-i T{ — , ■ ■ ■ ,-i 



dxi ' ' dxN 

For example, the ultrasecond-quantized identity operator corresponding to the 
identity operator in the bosonic case has the form 

OO OO oo oo „ „ 

Eb^Y^Y. fcuTi ■■■ dxi...dxkdyi... dy2M 

X &+(xi,si) • • ■b+{xk,Sk)B'^{yi,y2) ■ • • -B+(y2M-i, Z/2A/) 
X Symm^^_^^j^^ y^^^ (b~ixi, si) • • •S"(xfc, Sfc) 
X B^{yi,y2) ■ ■ ■ B^{y2M~i,y2M)) 

xexp^-^ J dxb+{x,s)b-{x,s)~JJ dydy' B+{y,y')B~{y,y') 

(14) 

which is the sum of the projections (jlip . Likewise, the ultrasecond-quantized 
identity operator in the fermionic case is given by 



oo oo 

jAsymm 

,M ' 

k=0 M=0 



and Symm in (I14p is replaced by Asymm. 
Consider the eigenvalue problem 



Hb,f^ = XEb^, £;$7^0, (15) 
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in the bosonic and fermionic cases. The following assertion was proved in [J]. On 

the subspaces and ^/f^™™ '^f space ,9' , the operators H b and Hp, 

respectively, acted upon by the projection onto the N -particle space (of the oper- 
ator (|2ip of the number of particles) coincide with the operator H^. Hence the 
eigenvalues A of problem (jisp in the bosonic and fermionic cases coincide with 
the corresponding eigenvalues of the operators p2)). If the commutators 
of the operators b~{x,s) and 6+(a;,s), as well as of B~'{x,y) and B^{x,y), are 
small of the order of then, according to [3], the asymptotics of solutions 
of problem (jlSp is determined by the symbol corresponding to problem (|15p . In 
the bosonic case, the pseudosymbol is defined in the form 

^Br{-)A-).B*{-),B{-)] 

E ■■■ dx,...dXkdy,...dy2M 

k,M=0 ' ' si=0 Sfc=0"^ 

X Si) • • • b*{xk,Sk)B*{yi,y2) ■ ■ ■ B* {y2M-i-,y2M)Hk+2M 

X Symm^^_^^y^....y^^,^ {b{xi, si) • • • b{xk, Sk) 

X B{yi,y2) ■ ■ ■ B{y2M-i,y2M))^ 

fc7j]j^ ■■■ dx[... dx'^, dy[... dy'^M' 

k',M'=o ■ ' s[=0 s'^,=a 

X b*{x[,s[) ■ ■■b*{x'^,,s'^,)B*{y[,y'^) ■ ■ ■ B*{y'^^,,^„y',^,,) 

X Symm^, ^, , y s'^) • • • 6(4,, 4,) 

xB{y[,y',)---B{y',^j,_„y',^j,))\. (16) 



The expression for the symbol in the fermionic case is similar except that Symm 
in (|16p is replaced by Asymm. The following identity holds for the pseudosym- 
bol (|16p in the bosonic case: 

^Br{-),b{-),B*i-),Bi-)]^^^^, (17) 

where H and ps are the second quantized operators, 

[ dx i,+ {x)( -^a)^-{x) + 1; II dxdy V{x,y)ij+{y)^+{x)^-{y)^-{x). 

(18) 



2ra ^ ' 2 
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Here pB depends on the functions h{x^ s) and B{y^ y')\ 

PB^Y^Y. fclM'ffc + 2Afv (^ dxdx' b{x,s)b*{x',s)^+{x)ii-{x') 

k=0 M=0 ■ ■ ^ ^ n J J 



^ M 

dyidy2 B{yi,y2)ip^ {yi)ip'^ {y2) 



dy[dy'2 B{y[,y'2)i^-{y[)i^-{y'2)^ """^{^ J ^^(^)V^" (^)) ' 

(19) 



where {x) and ip (x) are Wick ordered bosonic creation-annihilation opera- 
tors [7]. In the fermionic case, one has the similar identity 

where H and pp are the following second-quantized operators: 
H = [ dx ^+{x)( —^A]i,-{x) + 1 1 1 dxdyV{x,y)iA{x)iA{y)^j-{y)^-{x) 



= E E k\M\{k + 2M)\ [jj '^y^'^y^ B{y,,y2)^^^y,)i^+iy2)^ 



k=0 M=0 



OO OC ^ ^ 

X • • • / . . . / c?a;i dx'i . . . dxk dx'f. 

si=0 Sfc=0"' 

X b{xi,si)b*{x[,si) ■ ■ ■b{xk,Sk)b*{x'^.,Sk) 
X ij+{xi) ■ ■■i^+{xk)Poi^-{x'k) ■ • • V^-(x'i) 

X (^JJ dy[dy'2B{y[,y',)$-iy'^)$-iy[)^ . (20) 

Here '^"'"(j:) and '4'~{x) are the fermionic creation-annihilation operators, and Pg 
is the projection onto the vacuum vector of the fermionic Fock space. In general, 

the pseudosymbol of the ultrasccond-quantized operator L corresponding to an 
arbitrary second-quantized operator L is expressed by the formula [1] 

LsMb*{-)M-),B*{-),Bi.)]^^^pl^. 

In the space one introduces the ultrasecond-quantized particle number op- 
erators [3] 

OO OO OO OO 

^i^ = EE(^+2WfcrM'"> ^F = EE(^ + 2^^)nt"- (21) 

fc=o M=o k=a M=a 
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Accordingly, the pseudosymbol of the operator N b va. the bosonic case has the 
form 

^S = {EE ^^^J^E'-'E f ■■ [dx, ...dXk+2M 

'-fe=OM=0 '■ ■ si=0 Sfc=0-^ 

X Si) • • • b*{xk,Sk)B*{xk+l,Xk+2) ■ ■ ■ B*{Xk+2M-l,Xk+2M) 

X Symm^^ ^^^^^ {b{xi,si) ■ ■ ■ b{xk,Sk) 

X B{Xk+l,Xk+2) ■ ■ ■ B{Xk+2M-l,Xk+2M)) | 

E E fcTjj^E--- E ■■■ dz^...dZk'+2M' 

k'=OM'=0 ' ' s'=0 s',,=0 
1 fc' 

X b*{zi,s[) ■ ■ ■ b*{Zk>,Ski)B*{Zk' + l, Zk'+2) ■ ■ ■ B*{Zk'+2M'-l, Zk'+2M') 

X Symm^^...^^,^^^, {b{zi, s[) ■ ■ ■ b*{zk',s'k>) 

X B*{Zk'+l,Zk'+2) ■ ■ ■ B{Zk'+2M'-l,Zk'+2M'))^ 

(22) 

In the corresponding fermionic formula, Symm is replaced by Asymm. 

4 SYMBOL OF AN ULTRASECOND-QUANTIZED 
OPERATOR 

First of all, note that the above definition of pseudosymbol does not fully reflect 
the thermodynamic asymptotics, even though it complies with the Bogolyubov- 
Dirac rule saying that the creation-annihilation operators in the leading asymp- 
totic term should be replaced by c- numbers. In any case, we refer to the operator 
thus defined as the pseudosymbol. Let the operator H have the form 

L 



H = Y, I ■■■ I dx,...dxi i>+{x^)---i>+{xi) 
1=1 



1 1 

X Hi(^xi,...Ji;-i-^,...,-i-^Jij~{xi)---ij~{xi). (23) 

Then, in the case of ultrasecond quantization without creation-annihilation op- 
erators B^{x,y) for pairs of particles, the operators H and E defined above 
satisfy the identity 

H = EA, (24) 
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where A is an operator in ^ of the form 

L oo oc „ „ 

^ = ^ ^ • • • ^ . ■ ■ dxi . . .dxi 6+(a;i,si) • • •'0+(a;;,s;) 

1=1 si=0 s,=0"^ 

1 1 
/ 2 2 \ _ 

X H;LTi,...,a:;;-i^^,...,-i— 16 (xi,s;)---& (a;i). (25) 

If ultrasecond quantization also takes into account creation-annihilation op- 
erators for pairs of particles, then identity ([24]) remains valid, but the operator A 

has a more complicated form than (|25p . However, for pairs, the operator A has 
the form 



+ JJ dxdy B+{x,y)(^~^^{A, + Ay)^B-{x,y) 



B=0 ■ 



-j^ OO oo „ „ 
Sl=OS2=0-^-^ 

+ £ /// ^2;dyd^ + F(a;,z))6+(a;,s)i3+(y,z)B-(2/,z)6-(a;,s) 

+ dxdyV{x,y)B+{x,y)B-{x,y) 

+ l////d.*d.*„V'(.,„)B-(.,!,)g*(,.») 

X w)i?"(a;, z) + B^ {w,y)B^ {z,x) 

+ B^ {y, z)B^{w, x) + B^ {z, y)B^{x, w)^ . 

(26) 

If we replace the operators 13^{x,y) in (^5]) and ([^ by c-numbers, then 
we obtain the symbol corresponding to the asymptotics in the thermodynamic 
limit. 

5 BOSONIC CASE. BOGOLYUBOV FORMULA 

Consider a system of N identical bosons of mass m in a three-dimensional 
parallelepiped T with side lengths Li, L2, and L2- We assume that bosons 
interact with each other and the interaction potential has the form 

V{N'^'{x-y)), (27) 
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where V{^) is a compactly supported even function and x and y are the co- 
ordinates of the bosons on T. We assume the periodic boundary conditions 
along Li and impose the conditions that the derivatives are zero along L2. 
Note that the interaction potential (|?7|) depends on N in such a way that its 
range decreases with increasing particle number N , while the mean number of 
particles interacting with a given particle remains constant. 

An explicit expression for the ultrasecond-quantized operator H correspond- 
ing to the bosonic system in question under ultrasecond quantization with pairs 
is given earlier in this paper. As was discussed above, this ultrasecond-quantized 
operator satisfies the identity 

H = EA, (28) 

where E is the ultrasecond-quantized identity operator and A is an operator in 
the ultrasecond quantization space. It is easily seen that the operator 

A = JJ dxdy B+{x,y)(^-^{A, + Ay) + V{N'/''{x-y))^B-ix,y) 

+ ^JJJJ dxdydx'dy' V{N^'\x - y))B+ {x,y)B+ {x' ,y')B- (x,x')B- {y,y'), 

(29) 

where B^{x,y) and B~^{x,y) are the bosonic creation-annihilation operators 
for a pair of particles in the Fock space of ultrasecond quantization, satisfies 
identity (j28|) . By (f28|) . to find the asymptotics of the spectrum of the bosonic 
system in question in the limit as — > cxd, one needs to find the corresponding 
asymptotics for the operator (|29p. 

Since the function (l27l) multiplied by A^ weakly converges as A^ ^ cx) to the 
Dirac delta function, the second term in the operator ((29l) has the factor 1/N 
in this limit case. This means that, to find the asymptotics of eigenvalues and 
eigenfunctions of A, one can apply the semiclassical methods developed in [S]. 
The asymptotics of eigenvalues and eigenfunctions is determined by the symbol 
of the operator ; this symbol is called the true symbol of the ultrasecond- 
quantized problem. The true symbol corresponding to the operator ([29]) is given 
by the following functional defined for a pair of functions $+(a;, y) and $(a;, y): 

jr[$+(. ),$(•)]= JJ dxdy ^+ix,y)(^-^iA, + Ay)y{x,y) 

+ '^JJJJ dxdydx'dy' {NV{N^/^ix - y)))^+{x,y) 
X ^+{x',y')^{x,x')<^iy,y'). 

(30) 

Since the number of particles is preserved in the system for the functions $+(x, y) 
and ^{x,y), we obtain the condition 

dxdy ^+{x,y)<^>{x,y) = (31) 



12 



By the asymptotic methods in [5], to each solution of the system 

with condition (j3ip . there corresponds an asymptotic series of eigenfunctions 
and eigenvalues of the operator ([^5]) in the limit as ^ oo. It follows from the 
explicit form of the true symbol ([50]) that system ([5^ can be represented in the 
form 



f]$(x, y) = - — (A, + Aj,)$(x, y) 



+ jjdx' dy' {NV{N^^^{x - y)) + NV{N^^^x' - y'))) 
X ^+{x',y')<^{x,x')<^iy,y'), 



n^+{x, y) = -^{A, + Ay)<i>+{x, y) 



+ 2 jjdx' dy' {NV{N'^'^{x - x')) + NV{N^I^{y - y'))) 
X ^{x',v')^+{x,x')^+{y,y'). 

(33) 

Let Vq be the coefficient in the Fourier series expansion of the potential NV{N^/^x) 
on the parallelepiped (Li, ^2, -^2): 



e-w-'=Nv{\/N{x))dx, v^q^vq. (34) 



J2 JT 

The exact solution of system ([33|) is given by the functions 



= cos(fc2(x ~ y)), (35) 

^k,M = E ¥'fe.,ie^'^^"+^^e^'^--^\ (36) 

L1L2 



where fci = (7ii/Li,0, 0) and fc2 = (0, n2/£2, ?^3/-^2), with the eigenvalue 



n = —{kl + kl) + + V2k^, (37) 
m 



where the function ^pk2,i has the form 



b, 1 

2 
1 



1 

h^/m{P - kl) - (wo + V2k:) , , 
, bi — b^i. 

Vl-k2 + «/ + fe2 



(38) 



13 



(If vi 0, then ipk2J ^ 0-) 

The pair (fci, ^2) of vectors plays the role of parameters numbering various 
solutions of this system. The vector hki /m is equal to the bosonic system flow 
velocity along the capillary. The vector ^2 is the wave vector of the transverse 
mode. 

Note that ^ 00 as |Z| — > 00, since 

J^lJ^2 JNT ^1^2 JNT ^1^2 J 

(39) 

and consequently 

1 , . 

'Pk,i = -^, (40) 

whence it follows that the series ([55]) converges absolutely. 

We split the series (1551) into two parts, one with I < N^^^ and the other 
with I > N^/^. The first part of the sum converges as iV ^ oo modulo 0{N^^^^) 
to 

1 f7^ 7 dof 



' I /l2 1 aoi 



This readily follows from the change of variables v Nx = f in 

The second part of the sum converges to zero as 0{N^^/^) by Hence 
system ([33| supplemented with condition ((3T|) has the following family of solu- 
tions for /ci = in the limit as ^ oo: 



*fc (a;, y) = -rT2 coa{k{x - y)), 

^1^2 

y) = y~t2 Y1 '^f'''- exp(«?(a; - y)), 



(42) 



where k and Z are three-dimensional vectors of the form 

n2 ^3 ' 



\ ' L2' L2 

712 and 713 are integers, the ipk,i in (|28[) have the form 



"2..' = 2^ (l;<'^' - ' ± \/ - P) + v.) - V?) , (43) 

(the plus sign is taken for P > /c|, and the minus sign, for P < A:|), and finally Vq 
stands for the expression 

Vo = f dxV{x), (44) 
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the integral being taken over R^. The vector k in (|^^ plays the role of a 
parameter numbering various solutions of system ((33|) . (IST]) . The solutions (|42p 
are standing waves; there is no flow in the corresponding series. 

The leading asymptotic term of the eigenvalues for the series corresponding 
to the solution ([55)1 . is equal to N times the value of the symbol ([50]) on 
these functions: 



The asymptotics of eigenvalues and eigenfunctions (in particular, the terms 
following E}.^ }.^) is determined not only by system ([55)1 but also by the solu- 
tions of the variational system corresponding to the Hamiltonian system. The 
variational system for p3p has the form 



{n - X)F{x, y) = -|^(A. + Ay)F{x, y) 



+ jj dx' dy' {V{y/N{x - y)) + V{y/N{x' - y'))) 

X {G{x' ,y')^{x, x'My, y') + ^+ {x' , y')F[x, x')^{y, y') 
+ ^+{x',y')<^{x,x')F[y,y')), 



{n + \)G{x, y) - -:J^(A, + Ay)G{x, y) 
2m 



+ JJ dx' dy' {V{\/N{x-x')) + V{\/N{y-y'))) 

X {F{x', y')^+{x, x')^+(y, y') + $(x', y')G{x, x')^+{y, y') 
+ ^x',y')^+{x,x')G{y,y')). 

(46) 

To find the spectrum of quasiparticles, one should select solutions of the 
variational system satisfying the selection rule lOj for the complex germ for 
self-adjoint operators with real spectrum. 

Consider the case in which ^2 ^ 0. We substitute the solutions ([55)1 . (|36p 
into (|46|) and take into account symmetry. Then we see that the solutions of 
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the variational system have the form 

Gi{x, y) = ui^i(exp(i(-fci + k2)x + i{-ki + l)y) + exp(z(-fci + k2)y + i{-ki + l)x)) 
+ U2,i(exp(i(— fci — k2)x + i[—ki + 2k2 + l)y) 

+ exp{i{~ki -k2)y + i{-ki + 2fc2 + l)x)) , 
Fi{x, y) = -t;i,;(exp(i(A:i + k2)x + i{ki + l)y) + exp(z(/ci + k2)y + i{ki + l)x)) 
- V2J (exp(z(A;i - k2)x + i{ki + 2k2 + l)y) 

+ exp(i(fci - k2)y + i{ki + 2k2 + l)x)) 
+ X! wi,v{eyi-p{i{ki+k2 + l-l')x + i{ki+l')y) 

l'=^l,l+2k2 

+ cxp(i (fci + k2 + l- l')y + i{ki+ l')x)) , 

(47) 

where I ^ —k2 and the numerical coefRcients uu, U2.1, wi.i, V2^i, and w;^;/ are 
determined from an infinite system of equations. This system contains a closed 
subsystem of four equations for the coefficients Uij, U2j, fi,;, and V2^i, which 
can be rewritten in the form 

XX = MX, 



(48) 



where 



A = Ah ki{k2+l) 

m 



X 



and X is the column vector 

U2,l 
^1,1 
\v2aJ 

Equations (gT]) and ([48]) give the matrix M 

V2k2 + '"l + k2 



( 



M 



Bi + 



Vl + k2 + '"l-k2 



V2k2 + Vl + k2 



2{V0 + Vl^k2)V>k2,l 



{V2k2 + fi + fta) 
x{'Pk2,l + fk2,l + 2k2) 

{V2k2 + "i + fca) Of. I ^ 

„ Z r. \ ^V^Q + ^l + 3k2)'Pk2,l + 2k2 

\ XWk2,l + Vk2,l + 2k2) 



-B, 



Vl-k2 



V2k2 + Vl + k2 



where the Bi have the form 



'^l + k2 + ^l + 3k2 



V2k2 + Vl+k2 



-B, 



l + 2k2 



V2k2 



— {I - fca) + {vi-k2 + vi+k2)^k2,i Y 



Obviously, 



Vl + k2 + Vl + 'ik2 = Vl-k2 + Vl + k2 



(49) 
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uniformly in Z as oo (it suffices to make the change of variables (15^). 

Then the matrix M can be approximately represented as the block matrix 



M = 



C -ViE 
D ~C 



where E is the identity 2x2 matrix and 

Vl-k2 + Vl + k2 



Vi = 

We also introduce 



. _ Vl + k2 + V2k2 rr- _ Vl~k2 + VQ 

- 2 ' ' " 2 ■ 



The eigenvalue corresponding to Eq. have the form 

Afci,fe2,i = -2afci(fc2 + /) 



1/2 

(50) 



where 



a = , li = I + 2fc2- 

2m 



Before passing to the limit, one has 

1/2 



Xk.Md = -2afci' ± yiar + Vi - Vi+Y + ^ - ± 2(aZ" + Vi- Vi+)\V+\ 
Vi = Wi, = V; = — ^ — 

for k2 = 0. By formally setting k2 = 0, we arrive at Bogolyubov's well-known 
formula 

Here vi is the Fourier transform of the potential. 

(We assume that L2 is much larger than some standard length, say, the 
electron radius tq, and that although Li 3> L2, we can take a sufficiently large 
integer ni. In other words, L2/rQ — > 00 and Li/tq 00, but the vector ki = 
(ni/Li,0, 0) remains finite, since ni — > 00. 

In the language of nonstandard analysis, this means that L2 is an infinite 
(nonstandard) number, Li and ni are nonstandard numbers of higher order, 
and ki — (ni/Li,0,0) is a standard finite number. 

Then ^2 is equal to an infinitesimal nonstandard zero, fc2 = 0, and ki is a 
standard number.) 
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6 CASE OF A FERMI FLUID 

Consider the Hamiltonian system for fermions: 

\ 2m J 
+ '^N jj dx'dy' iV{x-y) + V{x' -y'))'S?+{x',y'Mx,x'W,y), 

m+ix,y)= f-^(A, + Ay))<^>+ix,y) 



2m 



2N II dx'dy' {V{x-x') + V{y-y')mx',y')^+{x,x')^+{y',y). 



(51) 

The functions y) and <i>(x,y) are antisymmetric and satisfy the normal- 

ization condition 

dxdy <^^{x,y)<^{x,y) ^ -. (52) 



(53) 



2 

Let us represent the interaction potential by a Fourier series: 

NV{x) = V Wpe*"^, Vp = — ^ / dx NV{x)e-'P'', Vp = 
p L1L2 J 

We seek the solution of system (|5ip , ([52]) in the form 

Kk^i^'V) = -l^e-*'=^(-+^) sin(fc2(^-y)), 

where fci, ^2, and / are three-dimensional vectors of the form 

\Li L2 L2 J 

and ni, n2, and are integers. The numbers ipk2.i should satisfy the condition 

After the substitution, we find that the eigenvalue is equal to 

V, = — {kl + kl) -f V2k2 - 
m 

and the ipk2.i have the form 

*&fc2,; , 1 n ;r~ u -kl) + {va-V2k2) 

2 2 V vi^k2 - vi+k2 
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Note that 
Set 



bk2,i — -bk2,-i- 
ibkid , 1 vi-k2 - vi+k2 



2 2 \vi^k2 - '"i+k2 
Then (pk2,i will be equal to —ipk2.-i- 



k2,l ■ 



^^(A.. + A,)F(a;,y) 



Consider the fermionic variational system 

+ '^^ JJ dx' dy' {V{\^{x - y)) + V{\^{x' - y'))) 

X {Gix\y'Mx, x'My', y) + ^+ [x' , y')F{x, x')^{y' , y) 
+ ^+{x',y')^x,x')F{y',y)), 

{n + X)G{x, y) - -^{K + ^v)G{x, y) 
2m 



-2N 



dx' dy' {V{\/^{x - x')) + V{VN{y - y'))) 



X {Fix', y')^+ix, x')^+iy', y) + ^{x' , y')G{x, x')^+{y', y) 
+ <^{x',y')<^+{x,x')G{y',y)). 

(54) 

Its solutions have the form 

Ql{x y) = Ui i(^Q^(-'''^+^'^)^+^i-^'^+^)y — g2(-'si+fe2)a+i(-fci+/)x-j 

_|_ ^^^i{-ki-k2)x+i{-ki+2k2+l)y _ ^i(-ki-k2)y+i{-ki+2k2+l)x-j 

Pl(^X, y) = Vi i{^e^^^^+*'^y-'+'^(*'^+^)y — ^i{kl+k2)y+i{ki+l)x^ 

_j_ ^2 ^(^Qi(ki-k2)x+i(ki+2k2+l)y _ ^i(ki-k2)y+i{ki+2k2+l)x-^ 
_|_ ^ y^^^i(^^i(ki+k2+l-l')x+i(ki+l')y ^i(ki+k2+l-l')y+i(ki+l')x-^ 
l'^l,l+2k2 

(55) 

where I ^ —k2 and the numerical coefficients ui^i, U2,i, wi,;, 'y2,i, and w;^;/ are 
determined from an infinite system of equations. This system contains a closed 
subsystem of four equations for the coefficients ui^i, U2^i, ui,;, and W2,z, which 
can be rewritten in the standard form 
The matrix M is given by 



M 



Bi + 



Vl-k2 



'Vl + k2 ~ V2k2 



Vl+k2 — V2k2 



Bl+2k2 + 



'"l + 3k2 



Vl-k2 ~ ^l + k2 

2 




2i{yi^k2 - Vo)<Pk2,i 
i{v2k2 - ""i+fta) 

\ X{'Pk2,l + 2k2 - 'Pk2.l) 



X(((9fc2,; + 2fe2 - 'Pk2,l) 
2i{vo — Vl+3k2)'Pk2,l+2k2 



VI -k2 







Vl+3k2 — 'Vl + k2 



Vl + k2 ^ 1^2*2 



Vl + k2 - V2k2 



-B, 



vi+ak2 

i+2fc2 ^ — J 
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and the Bi have the form 



Bi^ —{I - k^) + i{vi+k2 - vi-k2)^k2,i 
Zm 



V2k2 

2 



For ^2 = 0, the eigenvalues have the form 

A2J 



ft2 ;^2;2 
-^l.i = 'fcl + -7, , 

m 2m 



— Iki + 
m 



2m 



+ V1-V0 



This means that we choose the vectors ^2 = (0, 1/L2,0) andfci = (n/Li,0, 0). 
If n ^ 1 and hence ki ^ k2, then the fluid vefocity vector is mainly directed 
along the tube. As before, let Li = 00, L2 = 00, and n = 00 be nonstan- 
dard numbers, let Li ^ L2, and let ki be a standard number. Then k2 = 0, 
and (fk2,i takes nonstandard values for k2 = 0. However, if < Z < \/^M, 
where M < 00 (i.e., < Z < 00), then ipk2,i can be assumed to be a stan- 
dard number. For k2 = 0, taking into account the selection rule, we obtain the 
formula 

^2,1 = Iki 



2m 



+ V1-V0 



by analogy with the bosonic case and formula ([T]) . 

As was mentioned in Sec. [21 the behavior of the expression under the modulus 
sign is similar to that of the Landau curve for bosons. The criterion for ki has 
the form 



\k^ 



m 

< -TT mm 

h? I 



Vl - Vq 



2m 



similar to the Landau criterion (the vapor-destroying velocity). We set the 
parameter k2 to be an infinitesimal nonstandard number. However, it is not 
exactly zero, and the presence of a nonzero k2 results in a spectral gap. To 
compute the gap in the first approximation, one has to find the spectrum of the 
matrix M modulo 0(fc|). Then, by using the selection rules, one can determine 
whether there is a gap in the spectrum. 
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